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Abstract. We report the first results from the BOOMERanG experiment, which 
mapped at 90, 150, 240 and 410 GHz a wide (3%) region of the microwave sky with 
minimal local contamination. From the data of the best 150 GHz detector we find 
evidence for a well defined peak in the power spectrum of temperature fluctuations 
of the Cosmic Microwave Background, localized at £ = 197 ± 6, with an amplitude 
of (68 ± 8)/iKcmb- The location, width and amplitude of the peak is suggestive of 
acoustic oscillations in the primeval plasma. In the framework of inflationary adiabatic 
cosmological models the measured spectrum allows a Bayesian estimate of the curvature 
of the Universe and of other cosmological parameters. With reasonable priors we find 
tt = (1.07 ± 0.06) and n s = (1.00 ± 0.08) (68%C.L.) in excellent agreement with the 
expectations from the simplest inflationary theories. We also discuss the limits on the 
density of baryons, of cold dark matter and on the cosmological constant. 



INTRODUCTION 



Acoustic oscillations in the primeval plasma before recombination (z ~ 1100) 
produce a delicate pattern of horizon and sub-horizon size structures in the Cosmic 
Microwave Background. In the framework of inflationary adiabatic perturbations 
the angular power spectrum of these structures features a harmonic series of peaks, 
with the first peak at multipole I ~ 200. 

Here we present a high confidence measurement of the first peak resulting from 
the analysis of the multiband, high quality image of about 3% of the sky obtained 
in the long duration flight of BOOMERanG [1,2]. We show that the BOOMERanG 
image is a faithful representation of the CMB at angular scales smaller than 10°, 
and we shortly discuss the methods we used to make sure that instrumental and 
astrophysical contaminations are negligible in our results. We present the measured 
power spectrum, and we discuss the cosmo logical implications of this measurement. 

THE INSTRUMENT AND THE DATA SET 

The BOOMERanG payload and the important subsystems have been described 
in [5-9]. It is a scanning telescope measuring simultaneously eight pixels in the 
sky. Four pixels feature multiband photometers (150, 240 and 410 GHz), two pix- 
els have single-mode, diffraction limited detectors at 150 GHz and two pixels have 
single-mode, diffraction limited detectors at 90 GHz. The NEP of these detectors 
is below 200 \iKqmb\[s at 90, 150, and 240 GHz, and the angular resolution ranges 
from 10 to 18 arcmin FWHM. The continuous scan of the sky (±30° in azimuth 
at l°/s or 2°/s ) allows a wide sky coverage. Moreover, it encodes different spher- 
ical harmonic components of the sky brightness into different frequencies within a 
detection bandwidth (~ 0.2-20 Hz) free from 1/f noise and drifts. The instrument 
is flown aboard a stratospheric balloon at 38 km of altitude, to avoid the bulk of 
atmospheric emission and noise. Using a long duration balloon flight carried out by 
NASA-NSBF around Antarctica we mapped ~ 1800 square degrees in the region of 
the sky where Galactic contamination is minimal [10]. The instrument collected 57 
million samples for each detector during the scans. A pointing solution matrix (ac- 
curate to 3' rms) has been constructed using the data from the attitude sensors (sun 
sensors, laser gyroscopes, differential GPS). The beam response has been calibrated 
at ground, using a tethered thermal source in the far field of the telescope. This 
calibration has been checked in flight observing compact HII regions. The respon- 
sivity has been calibrated against the CMB dipole, which is visible along the scans 
as a ~ 3 mK linear drift. Systematic effects limit the precision of this calibration 
to 10%. For further analysis the data have been high-passed in the time domain, 
so that 1/f noise and drifts are reduced. In the same process, all the structures in 
the sky larger than 10° are effectively removed. A time-time correlation matrix has 



been constructed using an iterative method which effectively separates signal and 
noise in the data [11]. The maximum likelihood map and its covariance have been 
obtained from the data vector, the time-time correlation matrix and the pointing 
matrix as described in [13]. We use the HEALPIX [14] pixelization with 14' pixel 
side as a compromise between coverage of high multipoles and computation speed. 

THE MAPS AND THE POWER SPECTRA 

The maps in the four frequency channels have been published in [1]. In figure 
1 we increase the contrast and show the central part of a map obtained using a 
naive coadding in pixels of data from three of the 150 GHz channels. The 90, 
150, 240 GHz maps all present the same degree-size structures, with amplitude 
ratios characteristic of CMB anisotropics (see figure 2), thus demonstrating that 
the bulk of the anisotropy we measure is cosmological in origin. Dust contamination 
is evident only at lower Galactic latitudes, and mainly in the 240 GHz channel. 
At high Galactic latitudes the 410 GHz channel is a good monitor of low level 
emission from interstellar dust. We find a good correlation between this map and 
the IRAS/DIRBE map extrapolated at 410 GHz (model 8 in ref. [10]) and filtered 
in the time domain in the same way as our data. This correlation allows us to 
estimate low level contributions from IRAS-correlated dust to the CMB map shown 
in fig. 1. We find (see [1,15] for details) that this dust can produce at most 2% 
of the mean square anisotropy plotted in the 150 GHz map. We also studied 
the effect of point sources in the observed region. We used the PMN survey and 
the WOMBAT [12] software to extrapolate the radio sources flux to 150 GHz. 
The rms brightness fluctuation due to all the sources in our field corresponds 
to ~ 7fiK C MB, and has to be compared to the total measured sky fluctuation 
~ 80fiK C MB- Its contribution to the angular power spectrum has been computed 
and is reported in fig. 3. The maximum likelihood power spectrum of the central 
part of the map (1% of the sky) is computed as described in [13,1] and is plotted 
in fig. 3. A distinct peak at multipole £ ~ 200 is evident, and significant power 
is detected at all the multipoles up to the highest observed (£ = 625). In figure 
4 we plot the rms signal in the center part of the maps at 90, 150, 240 and 410 
GHz, computed as AT rms = \Jj2e(2£ + 1)q/4?t. The AT rms have been converted 
into brightness fluctuations using the measured spectral efficiency of the different 
freqeuency channels. The rms of the 410 GHz brightness has been estimated 
from the component correlated with the IRAS/DIRBE maps; the amplitude of 
the thermal dust spectrum has been adjusted to fit such rms. The amplitudes 
of the other theoretical spectra for Galactic emission have been scaled from [16] 
taking into account the different window functions of BOOMERanG and COBE 
and assuming an £~ 3 power spectrum. Once again it is evident that the data from 
the three channels more sensitive to the CMB are fit very well by the derivative 
of a 2.726K blackbody, while all the other reasonable spectra for millimeter wave 
emission originated in our galaxy cannot fit our measurements. 



COSMOLOGICAL INTERPRETATION OF THE DATA 



Having been able to exclude all the suspect sources of local contamination in 
the data of the power spectrum shown in figure 3, we can reasonably assume in 
the following that the detected power is in the CMB and comes from structures at 
the recombination epoch (z ~ 1100). In the framework of adiabatic inflationary 
models, acoustic oscillations in the primeval plasma at horizon and sub-horizon 
scales produce a distinctive series of peaks in the power spectrum of the CMB. 
The BOOMERanG data support this picture. The position, amplitude and width 
of the peak evident in fig. 3 are consistent with the general adiabatic inflationary 
scenario [17], while the simplest models based on topological defects do not fit the 
data as well [18]. Using a quadratic fit we find that the peak is located at multipole 
I = (197 ± 6) (lcr). This result can shift a little bit if one assumes skewed power 
spectra similar to the adiabatic scenario ones [19]. This position is consistent with 
a flat geometry of space, but it is not univocally related to the density parameter 
Q if we allow for a non vanishing cosmological constant [20]. The continuous line 
in fig. 3, a very good fit to the measured data, has parameters (fl, fib, fl\, h, n s ) 
= (1.035, 0.06, 0.425, 0.75, 1.0) where fib, fl\, fl are the density parameters for 
baryons, cosmological constant, and for the total mass-energy density respectively; 
h is the dimensionless Hubble constant, n s is the spectral index of the spectrum 
of primordial density perturbations. A full Bayesian analysis has been carried out 
[2] in order to measure the cosmological parameters and compute significant con- 
fidence intervals. In order to do this, we need to specify the prior distribution 
assumed for each of the parameters. This is especially important in the case of 
CMB power spectrum measurements, since an important geometrical degeneracy 
is present [21], so different combinations of the parameters produce very similar 
power spectra. Given a flat model with fl\ ~ 0.7, a closed model producing the 
same power spectrum can be found decreasing fl\ and h and increasing fib, and 
an open model can be found doing the reverse. The marginalized likelihood curve 
for fl will then depend on the density of models found in each direction along that 
particular path. Our 95% confidence intervals for fl range from (0.88 — 1.12) to 
(0.97 — 1.35) depending on the assumed priors and parametrizations [1,2]. This 
strongly suggests a flat geometry of the Universe, and at least implies, with 95% 
confidence, a curvature length larger than 2 to 3 times the Hubble length today. 
This measurement has very important cosmological consequences. Several indepen- 
dent measurements point to a matter density parameter significantly below unity 
[22-26]. As a consequence a different form of energy is needed to fill the gap be- 
tween the total energy density of the Universe and the energy density in matter. 
A cosmological constant is a good candidate, and has the appealing feature of ex- 
plaining the recent evidence coming from the observation of distant supernovae 
[27,28] for an acceleration in the expansion rate of the Universe . The fine tuning 
problem characteristic of the cosmological constant [29] can be overcome invoking 
quintessence [30-34] . In the same adiabatic perturbations framework, we constrain 
n s in the range (0.71 — 1.01) to (0.86 — 1.24) (95% confidence) again depending on 



the priors [2] . Our constraints on Q and n s are consistent and support the simplest 
inflationary scenarios [35-37]. The third parameter constrained by BOOMERanG 
is the density of baryons. In the power spectrum this parameter controls the rela- 
tive amplitude of the first peak to the second one (see e.g. [17]). From the power 
spectrum of fig.4, this ratio is larger than expected for a standard primordial nu- 
cleosynthesis (Vt b h 2 = (0.019 ± 0.002), la [38]), thus suggesting a physical density 
of baryons somewhat higher. Depending on the assumed priors, we obtain 95% 
intervals for Q b h 2 ranging from (0.019 - 0.045) to (0.026 - 0.048) [2]. This should 
be considered as a spectacular agreement, since we are comparing the density of 
baryons inferred from the physics of acoustic oscillations ~ 300000 years after the 
big bang to the density of baryons inferred from the physics of nuclear reactions a 
few minutes after the big bang. Moreover, the measurement techniques are com- 
pletely orthogonal and subject to completely different systematic effects. Should 
more precise measurements imply a more significant disagreement, new physics will 
be needed, and several hypothesis have been proposed already [39-45]. One way 
to improve the Q b consistency is to force the n s value to be smaller. n s ~ 0.9 is 
not the best marginalized value for n s , but the fit to the data is very good and 
still consistent with inflation (see e.g. [46]). In fig. 5 we plot the BOOMERanG 
constraints in the (Q m ,Q\) plane. Here Q m = Q b + &c, with Q c density parameter 
for non relativistic dark matter. It is evident that the data are fully consistent with 
a flat geometry of the Universe, represented by the Q m + Q\ = 1 line; however, the 
CMB anisotropy data alone cannot distinguish between the different contributions 
to the total mass-energy density. This degeneracy is efficiently removed once addi- 
tional data are considered, like the data from distant supernovae or the data from 
the large scale distribution of galaxies. The combination of our data with any of 
the two is quite powerful, and suggests a 'concordance' model with ~ 0.7 and 
Q m ~ 0.3 (see [2] for details). In ref. [47] all the available CMB power spectrum 
data and the power spectrum of the distribution of IRAS galaxies from the redshift 
survey data are fit simultaneously, confirming the results above. 

CONCLUSIONS 

BOOMERanG has produced a faithful, wide map of CMB anisotropy at angu- 
lar scales smaller than 10°. The power spectrum derived from the central part of 
the map (~ 1% of the sky) is consistent with the inflationary cosmological model 
with gaussian adiabatic density fluctuations. A careful analysis of the power spec- 
trum data, rigorously taking into account the priors and the degeneracies related 
to these parameters, has been carried out. We conclude that the geometry of our 
universe is close to flat, and that the primordial density flucutations have a nearly 
Harrison-Zeldovich power spectrum. The density of baryons inferred from these 
measurements is slightly higher than the standard nucleosynthesis value. Combin- 
ing the BOOMERanG data with data from the observations of distant supernovae 
or from large scale structure studies we obtain significant determinations of both 



tt m = (0.48 ± 0.13) and fi A = (0.66 ± 0.07) (la). All these results come from a 
single detector, and use the central part of the observed region. We are currently 
working on the combined analysis of all the 12 detectors sensitive to CMB. 
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FIGURE 1. The central part of the region observed by BOOMERanG at 150 GHz. The data 
have been high-pass filtered in the time domain to remove instrumental drifts. This process 
also removes structures in the sky larger than 10°. No further filtering is applied. The map 
uses HEALPIX pixelization with 7' pixel side. Spectral and amplitude arguments show that the 
structures present in the map are degree-scale anisotropics of the Cosmic Microwave Background. 



FIGURE 2. 3D and 2D scatter plots of the brightness measured by BOOMERanG at 90, 150 
and 240 GHz in each 7' pixel at high Galactic latitudes (b < -20°; 84205 data). The units are 
toKcmb, so that pure CMB fluctuations produce 2D scatter plots with slope 1. 
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FIGURE 3. Power spectrum of 1% of the sky from one 150 GHz channel of BOOMERanG 
(filled circles). The continuous line is the best fit adiabatic inflationary model. The open triangles 
represent upper limits to systematic effects in the data; the filled triangles are our estimate of 
the power from point sources in the observed region; the open circles represent our estimate of 
the power from IRAS-correlated galactic dust emission. Note that all these fluctuations must be 
added in quadrature. 




FIGURE 4. rms brightness fluctuations measured by BOOMERanG and spectra for CMB and 
local sources. The CMB curve has been computed from the derivative of a 2.726K blackbody. 
The thermal dust curve has been normalized to the IRAS-correlated component detected at 410 
GHz in the BOOMERanG data. 




FIGURE 5. 95% confidence level constraints on cosmological parameters in the (f2 m , Qa) plane, 
from the BOOMERanG power spectrum (continuous contour) and from the high redshift super- 
novae (dashed contour). 



